Численный анализ: Основы численного дифференцирования

Численное дифференцирование означает высокорисковый переход от бесконечной гладкости исчисления к дискретному, конечному миру цифровых вычислений. Мы заменяем бесконечно малое предельное значение на измеримый шаг $h$. Хотя теоретическая производная функции $f$ в точке $x_0$ определяется как $$f'(x_0) = \lim_{h \to 0} \frac{f(x_0 + h) - f(x_0)}{h}$$, компьютерные системы не могут напрямую вычислять предел. Вместо этого мы используем формулы конечных разностей, что приводит к количественно измеряемому штрафу, известному как ошибка округления.

1. Геометрия производной

Чтобы приблизительно вычислить $f'(x_0)$, мы смотрим на соседние точки. В зависимости от выбранного направления, мы получаем два основных выражения:

Формула прямой разности: Используется, если $h > 0$. Она смотрит вперед к точке $x_0 + h$.
Формула обратной разности: Используется, если $h < 0$. Она смотрит назад к точке $x_0 + h$ (где $h$ — отрицательное число).

В реальной инженерии, например, при вычислении длины дуги криволинейной траектории, мы часто полагаемся на эти приближения: $$L = \int_{0}^{48} \sqrt{1 + (f'(x))^2} dx = \int_{0}^{48} \sqrt{1 + (\cos x)^2} dx$$ Если функция $f(x)$ известна только в отдельных точках датчиков, численное дифференцирование является единственным путём вперёд.

2. Математическое вывод через интерполяцию

Чтобы приблизительно вычислить $f'(x_0)$, предположим сначала, что $x_0 \in (a, b)$, где $f \in C^2[a, b]$, и что $x_1 = x_0 + h$. Мы строим первый полином Лагранжа $P_{0,1}(x)$, определённый точками $x_0$ и $x_1$:

Шаг 1: Построение интерполирующей функции

$f(x) = P_{0,1}(x) + \frac{(x - x_0)(x - x_1)}{2!} f''(\xi(x))$

Шаг 2: Дифференцирование

Дифференцируя обе части уравнения и подставляя $x = x_0$, получаем фундаментальное соотношение: $$f'(x_0) = \frac{f(x_0 + h) - f(x_0)}{h} - \frac{h}{2} f''(\xi)$$

3. Член ошибки и сходимость

Член $-\frac{h}{2} f''(\xi)$ представляет собой нашу ошибку округления. Эта формула показывает, что точность составляет $O(h)$, то есть при уменьшении шага $h$ вдвое погрешность примерно также уменьшается. Однако следует быть осторожными: хотя меньшее $h$ уменьшает ошибку округления, в конечном счёте это увеличивает ошибку округления из-за вычитания почти одинаковых чисел в числителе.

🎯 Основной принцип: Конечная разность

Численное дифференцирование заменяет предел конечным отрезком. Качество нашего приближения зависит строго от размера шага $h$ и гладкости (второй производной) функции.

$f'(x_0) \approx \frac{f(x_0+h) - f(x_0)}{h}$ с границей ошибки $\frac{h}{2} \max|f''(\xi)|$

ВОПРОС 1

Какое математическое построение используется для вывода члена ошибки формулы прямой разности?

Первый полином Лагранжа интерполяции

Суммы Римана

Теорема о двух милиционерах

Ряды Фурье

ВОПРОС 2

Каков порядок точности формулы прямой разности $f'(x_0) = [f(x_0+h)-f(x_0)]/h$?

$O(h)$

O(h²)

$O(1/h)$

O(h³)

ВОПРОС 3

Рассмотрим функцию $f(x) = \ln x$ в точке $x_0 = 1.8$. Если $h = 0.1$, какова максимальная граница ошибки округления, если $|f''(x)| \leq 0.3086$ на интервале?

0.01543

0.03086

0.15430

0.00154

ВОПРОС 4

В контексте численного дифференцирования, почему очень маленький шаг $h$ (например, $10^{-16}$) может быть проблематичным?

Это приводит к значительной ошибке округления из-за катастрофического вычитания

Это увеличивает ошибку округления

Это делает функцию недифференцируемой

Полином Лагранжа становится неопределённым

ВОПРОС 5

Как называется формула $f'(x_0) = [f(x_0) - f(x_0 - h)]/h$?

Формула обратной разности

Формула прямой разности

Формула трёхточечного среднего

Правило трапеций

Вызов: Мастер-класс по дифференцированию и интегрированию

Синтез по модулю 4

Этот вызов проверяет вашу способность применять формулы дифференцирования, экстраполировать для повышения точности и интегрировать по сложным областям, как обсуждалось в Примере 1 и основных алгоритмах урока.

Вопрос 1

Используйте формулу прямой разности для приближённого вычисления производной функции $f(x) = \ln x$ в точке $x_0 = 1.8$ при $h = 0.1$, $h = 0.05$ и $h = 0.01$. Определите границы ошибки для каждого случая.

Решение:
Для $f(x)=\ln x$, $f'(x)=1/x$ и $f''(x)=-1/x^2$.
1. При $h=0.1$: $(\ln 1.9 - \ln 1.8)/0.1 \approx 0.54067$. Граница ошибки: $\frac{0.1}{2}|-1/1.8^2| \approx 0.01543$.
2. При $h=0.05$: $(\ln 1.85 - \ln 1.8)/0.05 \approx 0.54798$. Граница ошибки: $\frac{0.05}{2}|-1/1.8^2| \approx 0.00772$.
3. При $h=0.01$: $(\ln 1.81 - \ln 1.8)/0.01 \approx 0.55402$. Граница ошибки: $\frac{0.01}{2}|-1/1.8^2| \approx 0.00154$.

Вопрос 2

Используйте составное правило трапеций с $n = 4$ для приближённого вычисления $\int_{1}^{2} x \ln x \, dx$.

Решение:
$h = (2-1)/4 = 0.25$. Узлы: $x_0=1, x_1=1.25, x_2=1.5, x_3=1.75, x_4=2$.
Формула: $\frac{h}{2}[f(x_0) + 2(f(x_1)+f(x_2)+f(x_3)) + f(x_4)]$.
Значения: $f(1)=0, f(1.25) \approx 0.2789, f(1.5) \approx 0.6082, f(1.75) \approx 0.9793, f(2) \approx 1.3863$.
Расчёт: $0.125[0 + 2(0.2789+0.6082+0.9793) + 1.3863] \approx 0.6401$.

Вопрос 3

Примените экстраполяцию Ричардсона для определения $N_3(h)$ для $f(x) = \ln x, x_0 = 1.0, h = 0.4$.

Решение:
1. $N_1(0.4) = 0.841181$, $N_1(0.2) = 0.911608$, $N_1(0.1) = 0.953102$.
2. $N_2(0.4) = N_1(0.2) + \frac{N_1(0.2)-N_1(0.4)}{3} = 0.935084$.
3. $N_2(0.2) = N_1(0.1) + \frac{N_1(0.1)-N_1(0.2)}{3} = 0.966933$.
4. $N_3(0.4) = N_2(0.2) + \frac{N_2(0.2)-N_2(0.4)}{15} = 0.969056$.